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Abstract We review three Li problems. First, the Li 
problem in the Sun, for which some previous studies 
have argued that it may be Li-poor compared to other 
Suns. Second, we discuss the Li problem in planet host- 
ing stars, which are claimed to be Li-poor when com- 
pared to field stars. Third, we discuss the cosmological 
Li problem, i.e. the discrepancy between the Li abun- 
dance in metal-poor stars (Spite plateau stars) and the 
predictions from standard Big Bang Nucleosynthesis. 
In all three cases we find that the "problems" are nat- 
urally explained by non-standard mixing in stars. 
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1 The solar Li problem 

As illustrated in Fig. 1, the present day solar Li 
abundance (AiM = 1.05; I Asplund et alj|200a) is much 
lower than the mete oritic Li abundance (Au = 3.26; 
Asplund et al.ll2~009f) . This large depletion in the ob- 
served solar Li abundance by a factor of 160 rela- 
tive to the primordial solar system composition re- 
mains one of the most serious challenges of stan- 
dard solar models, which, as illustrated in Fig. 1, 
destroy only a minor amoun t (0.06 dex) of Li (e.g. 
D'Antona fc MazziteTlil Il984h . It is important to 
note that modern models using updated OPAL opaci- 
ties predict too much lit hium destruction during the 
pre-main sequence (e.g 



Piau fc Turck-Chieze 2002 



D'Antona fc MontalbaiJ 120031) . at variance with the 



observations of Li in solar-mass stars in young open 
clusters (see e.g. Fig. 6). Wh en the new low solar 
abundances ( Asplund et al.ll2009t ) are used the problem 
is reduced ( Sestito et al.ll2006l ) . but in general classical 
models have problems dealing with pre-main sequence 
convection, so other parameters are invoked to reduce 
the efficient l ithium destruction during the pre-main 
sequence (e.g. Ventura et al. 19981 ). 



1.1 Comparison to solar analogs 

A comparison between the Sun a nd solar analogs of 
one s olar mass and solar mctallicity ([Lambert fc Reddv 



2004J ) shows the Sun to be "lithium-poor" by a factor of 
10. This apparent pecu l iarity in the solar Li abundance, 
led lLambert fc Reddvi (|2004l ) to suggest that the Sun 
may be of dubious value for calibrating non-standard 
models of Li depletion. However, Pasquini et al. (1994) 
have shown that there are solar-type stars that have 



1 A X = \og{N x /N H ) + 12 



2 



a Li abundance as low as in the Sun. Nevertheless, 



Pasauini et al.l |l99J) comparison sample of solar-type 
stars span a wide range in stellar parameters (5400 K 

< Tefi < 6100 K, 3.6 < log g < 4.6, -1.6 < [Fc/H] 

< +0.2) and therefore they are not representative of 
one-solar-mass solar analogs. In Fig. 2, we restrict the 



comparison sample of iPasquini et al.l (j 19941 ) to only so- 
lar analogs within ±200K in solar effective tempera- 
ture, ±0.3 dex of the solar surface gravity and ±0.3 
dex of the solar metallicity. As we can see, these solar 
analogs seem to cluster in two groups, one with very 
high lithium abundances of Au ~ 2.4, i.e., 20 times 
higher than solar, and the other group with Li abun- 
dances as low as solar. Why are there no solar analogs 
with intermediate Li abundances? Why the number 
of solar analogs with low Li abundances is much lower 
than the number of analogs with high Li abundances? 
Could this be the reason whv lLambert fc Reddvl (|2004h 
only found stars with high Li abundances around one 
solar mass? The lack of stars with intermediate Li 
abundances in the Pasquini et al. (1994) solar analog 
sample, and the lack of stars with both intermediate 
and low Li abundance around one solar mass in the 



Lambert fc Reddvl (|2004r ) sample arc probably telling 
us that both samples have biases, perhaps due to a se- 
lection of stars mostly in one or two evolutionary stages, 
e.g., mainly young stars, which are known to have high 
Li abundances. 

The recent work by Pasquini et al. ( 2008| ) for solar 
analogs and solar twins in the solar-age open cluster 
M67, shows that solar twins (M67 stars around solar 
effective temperature) have Li abundance as low as so- 
lar, but stars 100 or 200 K hotter span a broad range 
in Li abundances. So, it is important that the tem- 
perature scale of the comparison sample is accurate; 
otherwise offsets of about 100 K may introduce a bias 
in the comparison between the Sun and stars. 

1.2 Comparison to solar twins 

Solar twins, stars with stellar parameters very similar 
to the Sun, are ideal targets to see if the Sun is normal 
(or not) in its Li abundance. Being so similar to the 
Sun, it is possible to obtain reliable stellar parameters, 
and provided the Sun and the twins are analyzed (and 
observed) in a consistent way, the temperature scale is 
accurate. Furthermore, being selected due to their sim- 
ilarity in colors and luminosity to the Sun, they should 
span a range of ages very close to solar, avoiding thus 
potential biases in the selection of comparison stars in 
only one evolutionary stage very different to the present 
Sun. 

Solar twins have been searched for a long time, and 
although interesting solar twin candidates like 16 Cyg B 



(HD 186427) were identified in the past, detailed anal- 
ysis showed that th ey were significantly dif ferent to the 
Sun (see review bv lCavrel de Strobell[l996h . 

When the first close solar twi n (18 Sco) was found 
( Porto de Mello fc da Siivalll997h . it seemed to have a 
Li abundance near s olar, but much better data (e.g. 



Melendez et al. 20061 ) showed that its Li abundance 
is actually three times higher than solar. One so- 
lar twin is certainly not an acceptable number for 
a comparison between the Sun and stars, so a large 
survey of solar twins was urgently needed. The two 
largest recent efforts for finding field solar twin stars 
are being undertak ing by the group of Y. Takeda (e.g 
Takeda et al.ll2007l) and by our group dMelendez et al 



200dlMelendez fc Ramirez l2007t lMelendez et al.ll2009at | 



Ramirez et al.l 120091 ). Importantly, whenever possible, 
we are obtaining very high S/N for our sample stars, 
because otherwise only upper limits can b e obtained for 
their L i abundances. Indeed, as shown bv lTakeda et al ' 



( 20071 ) in their Fig. 12, their data with S/N ~ 150 can 
only estimate upper limits for stars with Au < 1.5, i.e., 
they can only reliably determine Li abundances when 
they are three times higher than solar. 

Our solar twin survey has been performed mainly 
with the 2.7m telescope at McDonald observatory in 
the North and with the 6.5m Magellan Clay telescope 
at Las Campanas observatory in the South. We have 
also obtained some Keck+HIRES data in the North 
and VLT+UVES and HARPS data in the South. Our 
data has been taken at R = 60,000-110,000 and achiev- 
ing S/N = 200-1000. The first pilot data set taken at 
Keck resulted in the discovery of the second best solar 



twin, HD 98618 (jMelendez et alj|2006l) . about a decade 
after the discovery of the first solar twin 18 Sco. HD 
98618 seems to be a solar twin as good as 18 Sco, and, 
as this twin, it has also a Li abundance three times 
higher than solar. Learning from the experience of our 
pilot Keck observations, we improved our criteria to 
select the best solar twins, empirica lly adjusting our 
T c ff scale ( Ramirez fc Melendez! 2005 ) for an apparent 



zero-point problem (jCasagrande et al.l [2009) . This is 
probably the reason why our first solar twin run at 
McDonald was very successful. Besides confirming the 
solar twin nature of 18 Sco and HD 98618, we iden- 
tified t wo additional solar twins, H IP 56948 and HIP 
73815 (jMelendez fc Ramirez 1 l2007h . both with a low 



Li abundance similar to solar. HIP56948 remains to 
this date the star that most closely resembles the Sun, 
with a T e g similar to solar w i thin 1 K, as recently con- 
firmed bv lTakeda fc Taiitsul (|2009[) using Subaru+HDS 
observations. The year 2007 was very prolific for so- 
lar twin studies, bes ides the twins found by our group, 
Takeda et al.l (|2007l ) reported the discovery of the fifth 
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solar twin, HIP 110963. Interestingly, this solar twin 
has a high Li abundance of Au = 1.7. These five solar 
twins were starting to fill the Li desert seen in Fig. 2. 

The year 2009 has been even better, with many 
more solar twins found using our McDonald data 



(IRamirez et al.l 120091) and M agellan+MIKE observa- 
" 2009al) . which together have 



tions (Mclcndez et al 



found more than 30 stars very similar to the Sun. Thus, 
we are starting to have a large sample of solar twins 
for meaningful statistics, in particular for addressing 
the long-s tanding question on ch e mical peculiarities in 



the Sun ([M clcnd ez et al 
Gustafsson et al. 2009h . 



2009a; R amirez et al.l 2009 



In Fig. 3 we show the Li abundance for solar twins 
and solar analogs that have been analyzed for Li in 
our survey. Most stars from the Magellan run have al- 
ready been analyzed, but the McDonald Li analysis is 
just starting. Open circles show detections and filled 
circles upper limits. The Sun is also shown for com- 
parison. The first important point to note in this plot 
is that, unlike the lack of stars with both intermedi- 



ate and low Li in the lLambert fc Reddvl (|2004[ ) sample, 
and the lack of stars with intermediate Li abundances 



in the lPasauini et alj (|1994l ) sample, our sample nicely 
covers a broad range of Li abundances 0.6 < Au < 2.4, 
meaning probably that our sample is not affected by 
any significant selection bias. 

In Fig. 4 we show the Li abundances vs. T e ff of our 
solar twin and solar analog sample restricted to stars 
with mass within ±3% solar and [Fe/H] within ±0.1 dex 
solar. The Sun does not look peculiar on this plot. One 
star with relatively high Li abundance stands out. As 
shown in Fig. 5, where Li is plotted as a function of age, 
the high Li abundance of this star is due to its young 
age. This plot has a smaller number of stars than Figs. 
3-4 because here, besides the constraint to ±3% in mass 
and ±0.1 dex in [Fe/H], we show only stars for which 
we could determine reliable ages (> 2.5 sigma). This 
figure definitely shows that the solar Li abundance is 
not abnormal, at least not for a solar-metallicity solar- 
age one-solar-mass star, which at about 4.6 Gyr has 
already depleted a significant fraction of its original Li 
abundance. 

A comparison of our solar twin results to one-solar- 
mass stars in solar metallicity (±0.15 dex) open clusters 
(selected from the sample of ISestito fe Randichl liooil 
and including the results from Pasquini et al. 2008) is 
shown in Fig. 6. Again, the agreement is excellent, and 
reinforces a strong correlation between Li depletion and 
age for one-solar-mass stars 



No n-standard models (e.g.LM ontalb an fc Schatzman 



2000t ICharbonnel fc Talonl l2005l: IXiong fe Dene! 12009 
Do Nascimento et al.ll2009h can reproduce the observed 



data, as shown in Figs. 7 and 8. We are in the pro- 
cess of obtaining better Li abundances and ages for our 
sample stars, which can potentially constrain the range 
of initial rotational velocities of our solar twins. 

Based on the results shown above, we conclude that 
the solar Li abundance is not peculiar but a product 
of depletion due to non-standard mixing which affect 
both the Sun and the solar twins. 



2 Is Li depleted in stars with planets? 

Many works have compared stars with planets to field 
stars instead of stars without detected planets. This has 
be en usually done to have a large comparison sample. 
Gonzalez fc Laws! |2000l ) compared the Li abundance 



of stars with planets to the Li abundance of field stars 
with detectable Li, finding that (as previously believed 
for the Sun), stars with planets tend to have smaller Li 
abundances. This apparent peculiarity in the compari- 
son between stars with planets and field stars could be 
due to a selection bias, as we have shown in the previ- 
ous s ection for th e case of the Sun. Indeed, the same 
year, iRvanl (|2000l ) showed that stars with planets have 
Li abundances indistinguishable from field stars, when 
stars of the same temperature, age and composition 
w ere compared. 

Israelian et al. ( 2004h compared also planet hosting 



and field stars, showing that indeed they have the same 
Li abundance, except perhaps for stars around solar 
T e ff. It is important to note that this analysis is not 
fully consistent, as the comparison sample used by Is- 
raelian et al. was taken from the literature ( Chen et al 



20011 ). potentially having problems with different tem- 
perature scales. Also, as we have shown in the previous 
section, there are indeed many field stars around solar 
temperature showing a low Li abundance, thus, stars 
with planets probably do not have anomalous Li abun- 
da nces even around the solar T e g . 

Takeda fc Kawanomotol ( 2005 ) also compared stars 



with planets and field stars, showing that their Li abun- 
dances are identical at any temperature except for a 
small narrow range around 5850 K ± 50 K. This T e ff 
range is so narrow that the difference between stars 
with planets and field stars could probably vanish when 
more comparison st ars were included. Interestingly, 
Chen fc Zhaol (|2006l ) found that the discrepancy could 
occur at lower temperatures, around 5750 ± 50 K. 
Again, the lack of Li-poor solar analogs, could be mis- 
leading many authors to believe that stars with planets 
may be different to field stars around solar T c g, but our 
sample of solar twins with low Li abundances weakens 
those claims. 
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In the largest hom ogeneous Li study of stars with 
planets and field stars, Luck fe Heiter ( 2006| ) concluded 
that there is no discernible difference between planet 
hosts and comp arison star s. 

In Fig. 2 of Gonzalez ( 20081 ). planet hosting stars 
also seem identical to field stars, except in the nar- 
row T c ff range of 5850 ± 50 K. No conclusions can 
be made for stars with T c ff < 5800 K bec ause only 
three planet hosts are available. Note that the I Gonzalez 
( 20081 ) study is not homogeneous, as the Li abundances 
were obtained from different literature works with prob- 
ably different T e ff scales and different selection crite- 
ria. Thus, there could be potential selection biases, 
as we have shown for the case of the apparent low so- 
lar Li abundance. Indeed, we have just finished the 
analysis of 4 pla net hosting stars and 6 stars without 
detected planets (jMelendez et al.ll2009a), both a round 
the above temperature range where Gonzalez] ( 2008h 
found anomalously low Li abundances in planet hosts. 
As shown in Fig. 9, our planet hosting stars do not 
show anomalously low Li abundances with respect to 
our stars without detected planets. 

We conclude that stars with and without planets 
probably share similar Li abundances, at least to the 
level of precision of current analyses, although a defini- 
tive conclusion is not possible because Li depends on 
both mass and age (like we have shown for the Sun). It 
is important that in future studies the comparison be- 
tween planet hosts and stars without detected planets 
is performed for stars with similar parameters (mass, 
mctallicity, age), and both samples must be analyzed 
in a consistent way. Any claim for differences in Li be- 
tween stars with and without planets must be taken 
with skepticism unless reliable ages and masses are de- 
termined for planet hosts and the comparison sample. 



3 The cosmological Li discrepancy 

One of the most important discoveries in the study of 
the chemical composition of stars was made in 1982 
by Monique and Frangois Spite, who found an essen- 
tially constant Li abundance in warm metal-poor stars 
(Spite & Spite 1982), a result interpreted as a relic 
of primordial nucleosynthesis. Due to its cosmologi- 
cal significance, there have been ma ny studies devoted 



to Li in metal-poor field stars (e.g. iRvan et al 



Melendez fc Ramirezl2004lBoesgaard et al.l2005 



2005; lAsplund et al 



2007: Hosford et al 



Using the theory of big bang nucleosynthesis (BBN) 
and the baryon dens ity obtained from WMAP data 
(jDunklev et alj 120091 ) , a prim ordial Li abundanc e of 
A u = 2.7 2±° ° 5 R is predict ed (jCvburt et al.l l2008l see 
also e.g. Steigman ( 20071 )). which is a factor of 2- 
6 times higher than the Li abundance inferred from 
halo stars. There have been many theoretical stud- 
ies on non-standard BBN trying to explain the cos- 
mological Li dis crepancy by exploring the frontiers of 
vpl 

m 



^ — -ir^Ff — — 

new physics fe.g.jCoc et al.ll2009UJedamzik fc Pospelov 



Kohri fc Santosol l2009h . Alternatively, the Li 



1200 

problem could be explained by a reduction of the origi- 
nal Li stellar abundance due to internal processes (i.e. 
by stellar depletion). In particular, stellar models in- 
cluding atomic diffusion and mixing can deplete a sig- 



nificant fraction of the initial Li content (Richard et al 



l2005t Piau 2008h . although such models depend on 
largely unconstrained free parameters. Due to the un- 
certainties in the Li abundances and to the limited sam- 
ples available, only limited comparisons of models of 
Li depletion with stars in a broad range of mass and 
metallicities have been performed. 

In order to provide meaningful comparisons with 
stellar depletion models, precise Li abundances for a 
large sample of stars are needed. We have r ecently fin- 
ished such a study (jMelendez et al.l l2009bl ) . achieving 
errors in Li abundance lower than 0.035 dex, for a large 
sample of metal-poor stars (-3.5 < [Fe/H] < -1.0), for 
the first time with precisely determined masses in a rel- 
atively broad mass range (0.6-0.9 M Q ). 

Our recent work shows that Li is depleted in Spite 
plateau stars. As can be seen in Fig. 10, the spread of 
the Spite plateau at any mctallicity is much larger than 
the error bar. Also, there seems to be a correlation with 
T e ff at any probed metallicity. Actually the correlation 
is better when Li is plotted versus stellar mass (Fig. 
11), showing thus that Li has been depleted in Spite 
plateau stars at any mctallicity. In this figure we con 



front the stellar evolution predictions of iRichard et al 



( 20051) with our inferred stellar masses and Li abun- 
dances. The models include the effects of atomic dif- 
fusion, radiative acceleration and gravitational settling 
but moderated by a parametrized turbulent mixing; so 
far only the predictions for [Fe/H] =—2.3 are available 
for different turbulent mixing models. The agreement 
is very good when adopting a turbulent model of T6.25 
(see Richard et al. for the meaning of this notation) 
CharbonB iddfciiPiiiiitaaM Au = 2.64. The stellar Li abundances 



1999: 



2006l:IShi et al.l2007l: iBomfacio et al.^ use d above were obtained w ith the latest MARCS mod- 



20091 



Aokietal 



2009T) 



with ob- 
served Li ab undances at the lowe st [Fc/H] from as low 
A Li = 1.94 feoniiacio et"al 20071) to as high as A Li = 
2.37 (jMelendez fc Raimrezll2004l ) . 



els (jGustafsson et al.l 12008). but if we use instead the 
Kurucz convective overshooting models, then the re- 
quired initial abundance to explain our observational 
data would correspond to Au = 2.72. 
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Our results imply that the Li abundances observed in 
Li plateau stars have been depleted from their original 
values and therefore do not represent the primordial Li 
abundance. It appears that the observed Li abundances 
in metal-poor stars can be reasonably well reconciled 
with the pr edictions from stand ard Big Bang nucleosyn- 
thesis (e.g. Cvburt et al. 2008| ) by means of more real- 



istic stellar evolution models that inclu de Li depletion 



throu gh diffusion and turbulent mixing ([Richard et al 



20051) . We caution however, that, although encourag- 



ing, our results should not be viewed as proof of the cor- 
rectness of the diffusion+turbulence models of Richard 
et al. models until the free parameters required for 
the stellar modeling are better understood from basic 
physical principles. In this context, new physics should 
not be discarded yet as a solution of the cosmologi- 
cal Li discrepancy, as perhaps the low Li-7 abundances 
in metal-poor stars might be a signature of supcrsym- 
metric particles in the early universe, which could also 
explai n the Li-6 detections claimed for some metal-poo r 
stars (|Asplund et al.ll2006t IXsplund fc Melendezll2008| ). 
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Fig. 1 The s tanda rd solar model (e.g. 

iD'Antona fc Mazziteiul 119841 ) can not reproduce the 
large decrease in Li abundance at the present solar age. 
The Sun is represented by the symbol 
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Fig. 3 Magellan and McDonald solar analogs and twins 
covering a relatively broad range in stellar parameters (5640 
K< Tea < 5920 K), 4.31 < log g < 4.56, and -0.45 < [Fe/H] 
< +0.45). Open circles are detections and filled circles are 
upper limits 
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Fig. 2 Li in the solar analog sample of IPasauini et al.l 
l| 19941 ) restricted to stars with T e ff within 200 K solar, and 
both log g and [Fe/H] within 0.3 dex solar. Open symbols 
are detections, filled symbols are upper limits 
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Solar twins (M=1.00 = 0.03 M e , [Fe/H] = 0.0 = 0. 1 ) 
(Melendez et al. 2009, Ramirez et al. 2009) 
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Fig. 4 Magellan and McDonald solar twins with one solar 
mass (± 3%) and solar [Fc/H] (± 0.1 dcx). Symbols as in 
Fig. 3 
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Fig. 5 Li as a function of age for our Magellan and Mc- 
Donald solar twins with one solar mass (± 3%) and solar 
[Fe/H] (± 0.1 dex). Symbols as in Fig. 3. Only stars with 
reliable ages (> 2.5 a) are shown 
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Fig. 7 Comparison of non-standard solar models to field 
and open cluster stars 
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Fig. 6 Li for our solar twins with one solar mass (± 
3%) and solar [Fe/H] (± 0.1 dex), and for one-solar- 
mass stars in sol ar metallicity (±0 . 15 de x) open clus- 
ters selected from ISestito fc Rand ich (2005), although for 
M67 we used the sample of IPasauini et al l|200ct ). Field 
stars are shown as circles while open clusters (NGC 2264, 
IC2602/IC2391, Pleiades, Blanco 1, NGC6475, M34, Coma 
Berenices, Hyades, NGC762, M67) with triangles 
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Fig. 8 Non-standard models of Li depletion by 
ICharbonnel fc Talonl (|2005l ) for different initial rotation ve- 
locities 
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Fig. 9 Li in stars with (squares) and without (trian- 
gles) detected planets from the solar analog sample of 
IMelendez et al.l (|2009al ) . Open symbols are detections and 
filled symbols are upper limits 
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Fig. 10 Li abundances vs. T e fi for our sample of metal- 
poor stars in different metallicity ranges. The spread at any 
given meta llicity is much larger th an the error bar. Figure 
taken from IMelendez et alj lj2009bh 
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Fig. 11 Li abundances as a function of stellar mass in 
different metallicity ranges. Models at [Fe/H] = —2.3 in- 
cluding diffusion and T6.0 (short dashed lin e), T6.09 (dot- 
ted l ine) and T6.25 (solid line) turbulence (I Richard et al.l 
l2005h are shown. The models have been rescaled to an 
initial ALi=2.6 4 (long dashed line). Figure taken from 
IMelendez et al.l (|2009t 



